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Abstract 

Although the process of endocytosis of the low density lipoprotein (LDL) 
macromolecule and its receptor have been the subject of intense experimental 
research and modeling, there are still conflicting hypotheses and even conflicting 
data regarding the way receptors are transported to coated pits, the manner by 
which receptors are inserted before they aggregate in coated pits, and the display of 
receptors on the cell surface. At first it was considered that LDL receptors in human 
fibroblasts are inserted at random locations and then transported by diffusion 
toward coated pits. But experiments have not ruled out the possibility that the true 
rate of accumulation of LDL receptors in coated pits might be faster than predicted 
on the basis of pure diffusion and uniform reinsertion over the entire cell surface. It 
has been claimed that recycled LDL receptors are inserted preferentially in regions 
where coated pits form, with display occurring predominantly as groups of loosely 
associated units. Another mechanism that has been proposed by experimental cell 
biologists which might affect the accumulation of receptors in coated pits is a 
retrograde membrane flow. This is essentially linked to a polarized receptor insertion 
mode and also to the capping phenomenon, characterized by the formation of large 
patches of proteins that passively flow away from the regions of membrane 
exocytosis. In this contribution we calculate the mean travel time of LDL receptors to 
coated pits as determined by the ratio of flow strength to diffusion-coefficient, as 
well as by polarized-receptor insertion. We also project the resulting display of 
unbound receptors on the cell membrane. We found forms of polarized insertion 
that could potentially reduce the mean capture time of LDL receptors by coated pits 
which is controlled by diffusion and uniform insertion. Our results show that, in spite 
of its efficiency as a possible device for enhancement of the rate of receptor trapping, 
polarized insertion nevertheless fails to induce the formation of steady-state clusters 
of receptor on the cell membrane. Moreover, for appropriate values of the flow 
strength-diffusion ratio, the predicted steady-state distribution of receptors on the 
surface was found to be consistent with the phenomenon of capping. 

Keywords: LDL receptor dynamics, Polarized receptor reinsertion. Retrograde-flow 
model, Mean capture time. Surface-display capping 
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Background 

Endocytosis generally refers to the process by which substances are internalized 
through the cell membrane. Receptor-mediated endocytosis (RME) is a highly special- 
ized kind of endocytosis in which large protein molecules called receptors project from 
the cell membrane and couple selectively to ligands such as low density lipoproteins 
(LDL) [1]. Coated pits and their associated receptors have been studied most exten- 
sively in cells grown in culture, and the LDL receptor is the one on which a majority of 
experimental research has been conducted. The ligand-receptor complexes aggregate in 
specialized cell-membrane formations called coated pits [2-5]. When these close, they 
form vesicles which then transport the ligand-receptor complexes to the interior of the 
cell. Once the ligand-receptor bindings are separated, the ligands are degraded at the 
lysosomes. In some experimental systems, including the LDL system, the receptors are 
returned to the cell surface for further endocytic cycles, and the ligands are degraded at 
the lysosomes. Experimental results reveal that RME occurs in virtually all eukaryotic 
cells except the mature erythrocyte [6] and that it also provides an entrance mechanism 
for cell nutrients such as low density lipoproteins (LDL) [7], hormones and growth 
factors [8], serum transport proteins and antibodies [9], as well as toxins and lectins 
[10] and even viruses [11]. 

The low density lipoprotein (LDL) particles are cholesterol-transporting macromole- 
cules that are produced in the liver and circulate in the plasma. Membrane reconstitu- 
tion in human fibroblasts requires the assimilation of cholesterol. The efficiency of this 
process depends on the rate of endocytosis of the LDL macromolecule and its receptor. 
It is thought [12] that a severely depleted number of LDL receptors promotes high 
levels of circulating cholesterol because LDL internalization requires the receptor- 
binding stage. Deficiencies in the LDL cycle are known to be responsible for the 
ailment known as familial hypercholesterolemia, which promotes atherosclerosis and 
the incidence of strokes and coronary disease [13]; this being the reason why so much 
experimental research has been conducted on the internalization of the LDL receptor. 
A substantial portion of research aimed at the elucidation of the cholesterol-uptake 
process has addressed the characterization of the dynamics and display of the LDL 
receptor. This has produced a large pool of data and a sound conceptual framework 
which sustains the theoretical exploration of aspects of the dynamics mentioned above. 

It is generally thought that the maintenance of receptors on the cell surface is due 
primarily to receptor recycling rather that to de novo synthesis. These data come from 
experiments using cycloheximide to block protein synthesis. In this experimental setup, 
it is observed that the number of LDL receptors on the cell surface remain roughly 
constant for at least six hours [7]. Evidence of receptor internalization and reinsertion 
in unblocked systems would sustain the assumption that a steady-state concentration 
of receptors is maintained at the cell surface. The time receptors spend in the interior 
of the cell is negligible. The basis of this claim is the apparently undetectable pool of 
receptors inside the cell during endocytosis [14]. The transit time for an LDL receptor 
from binding on coated pits to reappearance in the membrane, and found it to be on 
the order of 15 seconds [14]. Based on these ideas, we will abide by the assumption that 
internalization and recycling of LDL receptors in human fibroblastic cells maintains the 
surface concentration of LDL at a steady state. Moreover, coated pits include 1% of the 
cell surface (coated pits include 2% of the cell surface at 4°C [15,16], but when the 
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temperature is raised to 37°C the number of coated pits on the surface is halved [17]). 
Anderson et al. [17] assert that coated pits tend to be linearly aligned over intracellular 
fibers. Hence, it can be assumed that the geometrical arrangement of coated pits on the 
cell surface of human fibroblasts can be reasonably approximated by means of a dilute 
and ordered system of sinks distributed on a two dimensional surface [18]. Guadorov 
et al. [19] reported that coated pits tend to form repeatedly at defined sites while other 
regions are excluded. Therefore, we can assume that coated pits are infinitely long- 
lived traps distributed in a dilute and ordered form over the cell surface. 

It is also known that on cultured human fibroblasts, receptors for certain ligands (e.g. 
insulin, epidermal growth factor and oc - 2-macroglobulin) cluster in coated pits only 
after exposure to the ligand [20], while receptors for LDL cluster in coated pits 
independently of ligand binding [7] . This feature of the LDL receptor pathway makes it 
a particularly attractive candidate for mathematical modeling, since initially we can 
ignore the details of the ligand-receptor interaction and still study the recycling of the 
receptor and its interaction with the coated pit. Preliminary mathematical models of 
RME are aimed at calculating the rate at which diffusing particles (receptors), which 
have been inserted uniformly all over a certain 2-dimensional space (the cell surface) 
hit traps (coated pits). This rate is known as the diffusion-limited forward-rate constant 
[21] and is denoted here by /crf+. This rate can be calculated as the flux of particles into 
a trap, divided by the mean particle concentration [22]. In the two-dimensional case, 
for a circular sink of radius a, kii+ is defined by means of the equation 

InaDdC . . , 

where C{r) is the steady-state radial distribution function of receptors not bound to 
coated pits, with D> 0 their diffusion coefficient and (Q the receptor concentration 
averaged over all the diffusion space [21], that is. 



(C) = ^ / 2nrC{r)dr. (2) 
nb 



The constant /cj+ times the number of traps per unit area gives the probability per 
unit time that a diffusing particle hits the trap. It is linked to the mean time for a 
particle to hit a trap (mean capture time) [21] by the relation, 

where p is the number of coated pits distributed per unit area on the cell surface. 

If LDL receptors are trapped in times less than, or comparable to the average lifetime 
of a coated pit (<5 min), the mean capture time is close to what it would be if the traps 
were infinitely long-lived. Moreover, since we are dealing with a dilute system of traps, 
we can follow the method of Berg and Purcell [23] and represent the real multiple-trap 
problem by means of single coated pit of radius a, surrounded by a reference annulus 
n of outer radius b, i.e. 

n = {{r,0)\O<a<2n,a<r<b}, (4) 



where b is calculated through 
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(5) 



Receptors are projected into the reference annulus fl according to an insertion rate 
function 5(r), and diffuse afterwards until they reach the boundary of the sink at r = a 
and are trapped. A steady-state concentration of diffusing particles will be maintained 
if the number of particles lost to the trap is the same as the number inserted per unit 
time. Then, if Cjs(r) denotes the steady-state receptor concentration at a distance 
r from the center of the sink we will have 

DV'^Cds{r) + S{r) = 0 for a<r<b, (6) 
so that the inner boundary at r = « absorbs, i.e. 

Cds{a) = 0 (7) 
and the outer boundary at r = reflects, i.e. 

9Qs 



dr 



0 at r = ^, (8) 



that is, there is no net flux of receptors across the outer boundary of fl. Solving for 
CdJj) we obtain 

r / uS{u)du 

Cdsir) = ^-^-^^ dz. (9) 

Denoting by means of CdJj) the form of obtained by setting S{r) = S, and with 

S being a constant, this means that whenever receptors are inserted uniformly over the 
entire reference annulus D we obtain 

Cduir)= In - (10) 

and if kdu + denotes the associated forward rate constant we have, 
_2jmDdCdu 

Then, for this case equations (1) through (3) yield 

'"'iDib^-a^)' 8D ■ ^^^^ 

Goldstein et al. [24] assert that whether traps are ordered or disordered makes only a 
small difference in the rate at which a trap captures receptors, given that the traps are 
distributed over the entire surface of the cell in a dilute way. Goldstein et al. [21,25] 
dealt with two conflicting hypotheses about the coated-pit recycling process. In either 
case, the coated pit effectively has a finite lifetime. However, they found that the model 
of equations (6)-(12) obtained under the assumption that sinks have infinite lifetimes 
give good approximations for kd+ in the experimental system for receptors of low 
density lipoproteins on human fibroblastic cells. Nevertheless, they also concluded that 
the experiments do not rule out the possibility that the true rate of accumulation of 
LDL receptors in coated pits could be faster than predicted on the basis of pure 
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diffusion and uniform reinsertion all over the cell surface. A mechanism that has been 
proposed by experimental cell biologists which might affect the accumulation of recep- 
tors in coated pits is convection. Evidence of a membrane flow comes from the obser- 
vation that, when the leading edge of a moving cell comes in contact with a small 
particle, the particle often adheres to the cell surface and is transported backward 
[26-28]. The motion of such adherent particles is consistent with particles undergoing 
Brownian motion in the presence of a constant membrane flow [29,30]. It has been 
hypothesized that this flow originates when membrane components like specific recep- 
tors and lipids, which had been internalized from the cell surface, are recycled to the 
leading edge and contribute to the extension of the cell boundary [29]. Considering that 
a rapid rate of membrane internalization is observed in endocytosis, a significant 
amount of recycled membrane should be made available for delivery at the leading 
edge. This delivery would be restricted to just those molecules that are internalized by 
coated pits [14,31,32], thus initiating the transport of membrane components away 
from the leading edge and towards the rear of the cell [26-31]. In addition to the 
specific delivery of receptors at the leading edge of the cell, their sorting within specific 
targets along the endocytic pathway have been proposed as components of a mechan- 
ism for the control of the internalization of specific ligands [18,33-35]. Generally, 
eukaryotic cells show the ability to efficiently control the traffic of many types of 
macromolecules, targeting proteins and lipids for a variety of destinations, based on 
their types and the current needs of the cell. This traffic control is crucial for the main- 
tenance of cellular structure and function [36,37]. This makes it reasonable to assume 
that reinsertion of LDL receptors could be arranged in such a way that their trapping 
rate might be adapted to suit specific metabolic requirements. Moreover, it has been 
claimed [38] that in giant HeLa cells, receptors are not inserted uniformly into the 
plasma membrane but rather at the periphery of the cell. Robenek and Hesz [39] 
claimed that their experiments with LDL particles bound to colloidal gold provided the 
first clear demonstration of the sequential clustering of receptors near coated pits in 
plaques or loosely associated groups, and not in the form of widely dispersed individual 
units. They concluded that this effect is produced when recycled LDL receptors are 
inserted preferentially in membrane regions where coated pits form. According to 
Wofsy et al. [40] plaques could be conceptualized as receptor clusters surrounding 
coated pits. The idea of a retrograde flow is essentially connected to the concept of a 
polarized reinsertion of receptors. However, neither the influence of this insertion mode 
on the aggregation rate of LDL receptors in coated pits, nor the resultant steady-state 
surface aggregation patterns have been formally studied. 



Methods 

For the aims of the present research we extend the formal framework presented in 
Echavarria-Heras et al. [41] for the purpose of including diffusion aided by a retrograde 
flow, plus polarized receptor insertion as a combination of mechanisms that influences 
both the mean capture time and the surface aggregation patterns of unbound LDL 
receptors. First, we conceptualize a partitioned receptor-insertion mode in which the 
cell membrane is divided into disjoint regions, each one receiving a fraction of the total 
number of recycled receptors. In the general form of this mode, different numbers of 
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receptors are projected over the different regions partitioning the cell surface. Polarized 
insertion, moreover, corresponds to the case in which receptors are recycled into two 
disjoint regions having abruptly contrasting rates. This partitioned mode can be used 
to obtain arrangements — which can be seen as special cases — comprising uniform [7], 
locally uniform, peripheral [38], and plaque forms [39,40], as well as a generalized radi- 
ally symmetric mode of the form described in [41]. 

In order to formalize the idea of a partitioned insertion mode, let's consider the refer- 
ence annulus D. = {{r,d)\ a< r <b, - n < 6 <n}, letting be the number of particles 
inserted in D, m and a being real numbers satisfying \<m< bla, 0 <a<Ji, we define 
the disjoint regions flc{m, a),flp{m, a) and 0.q{m., a) through 

D.c{m,a) = {{r,d)\ a < r < b, -n + a < 6 < n-a}, (13) 
np{m,a) = {{r, d)\ ma < r < b, n-a <d<n + a}, (14) 

Oqlm, a) = {{r, 6)\ a < r < ma, n-a <d<n+a}. (15) 
Then, we have 

D, = D,c{yn,a) yj D,p{m,a) yj D.q{m,a). (16) 

In order to characterize partitioned receptor insertion, we conceive an insertion rate 
function S'^^{c,p,q,m,a) that distributes receptors over the disjoint regions Dc(w, a) 



Figure 1 The geometry of the model, a) A circular trap of radius a (the coated pit) is encircled by an 
annulus of radius b (the reference region O associated with a coated pit). LDL receptors originally inserted 
at a point (f, 6) inside the reference annulus Q, move afterwards by convection and diffusion until they are 
trapped in coated pits, b) Receptor insertion occurs according to a partitioned insertion rate function 



a 



b 




S'^(c,p,q,w,a), which sorts receptors at distinct rates S'^(m,a), 
the disjoint regions Oc, Op and Qq. 



(m,a) and 5'^{m,a) linked respectively to 
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0.p{m, a), and Oqim, a) (Figure 1) through the respective non-negative and continuous 
functions S''f{m, a), S''p{m, a) and S''^{m, a). Formally, 

[S[\m,a) {r,d) eD,{m,a) 
S'-%c,p, q, m,a)=l Sfim, a) (r, 6) e a) , (17) 

with 

JL.Sf (m, a)dn + SSnXi^^ «)^^ + Il'r^,^^ (m, a)dn = In- (18) 

The formal properties of the partitioned insertion mode are elaborated in Appendix 
A. Table 1 summarizes the possible forms that S'^^(c,p,q,m,a) can achieve. Moreover, 
S'^^{c,p,q,m,a) can be handily characterized by using the functions Sc{m,a), 8p{m,a) 
and Sq{m, a), which respectively stand for the proportion of In that S''f{m, a), S'^^{m, a) 

and S'^{m,a) insert over the disjoint regions Dc(m,a), D.p{m,a) and D,q{m^a). For a 
and m fixed, equations (A4) through (A6) define the proportions Sc{m, a), Sp{m, a), and 
Sq{m, a). Whenever a pairwise combination of insertion proportions is vanishing, the 
associated form of the partitioned mode is polarized. Otherwise, the partitioned mode 
would be uniform all over D or might be locally uniform, thus sorting receptors at 
dissimilarly constant rates over Dc(w,a), D,p{m,a) and D,q{m,a). Although, for ease of 
presentation, we here partitioned fi into three disjoint regions, the form given by equa- 
tion (18) could be extended to include a greater number of disjoint regions forming £1. 

Let Cxmai^t S) denote the steady-state concentration of receptors in D.. In Appendix B 
we review the formal steps which establish that if receptors are inserted in D, in accord- 
ance with the insertion rate function S'^^{c,p,q,m,a) and afterwards move by diffusion 
aided by a retrograde membrane flow, then if A = v/2D, where v is the strength of the 
retrograde flow and D stands for the diffusion coefficient (cf. Eq. Bl), we will have 

CA^a(r,0) = GA„„(r,0)e^"°'', (19) 
where Gxmair, 6) is the solution to the boundary-value problem, 

^ (jXma[r,o)-A Gxma[r,0) = =; 6 , (20) 



Table 1 Different forms of the partitioned insertion rate function S'^Cc, p, p, m, a) : c, p 
and q stand for positive constants 



Mode notation 


Mode designation 


Mode defining conditions 




Eq. 


S\, 0, 0, bla, 0) 


c-uniform 


5?(^0)=c,S;«(^,0). 


= 0,5^(^,0). 


= 0 


A22 


S'®(0,c,c,m,7T) 


pq-uniform 


Sf{m,n) = 0,Sf{m,n) 


= S>,7,) = 


c 


A24 


S'^{c,c,c.m,d) 


cpq-uniform 


5f{m,a) = S|f(m,a) = 






A26 


S'^(c, p, q, m, a) 


cpq-locally uniform 


5f{m,a) = c,5f{m,a) 


= P,S^im,a) 


= q 


A28 


S%,p,q, m,n) 


pq-locally uniform 


Sf(m,n) = 0,Sf {m,n) 


= P,Sf(m,n) 




A30 


S%, 0, q, m,n) 


q-plaque form 


Sf{m,n) = 0,Sf {m,n) 


= 0,Sf(m,n) 


= q 


A32 


S'®(0,p, 0, m,n) 


p-peripherial 


5f{m,n) = 0,S^'*(m,7T) 


= P,Sf{.m,n) 


= 0 


A34 


S'®(0,p, 0, m, a) 


p-polarized 


Sf{m,a) = O.Sf{m,a) 


= p.S'^{m,a) 


= 0 


A36 


S'®(0,p, q, m,a) 


pq-polarized 


Sf{m,a) = 0,Sf(m,a) 


= p,S'^(m.a) 


= q 


A38 


S'^(0, 0, q, m, a) 


q-polarized 


5f{m,a) = 0,S|f(m,a) 


= 0,S^{m,a) 


= q 


A40 
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Ga«,„(«,0) = O, (21) 
Gxma{r.e) = Gxma{r,-e), (22) 
<Ji»,a('", ^) periodic in 0, (23) 



dr 



r = i-A COS e Gxma{b, 6) 6^'^''™'^ " dB = 0, . (24) 



The associated forward rate constant is denoted by means of k^a + and is given by 



nb J J S {c,p,q,m,a) rdrdd 

kxma+ = " ^f, l-jj , (25) 

C\ma{r, 0) rdrdd 



and the associated mean capture time Txma is given by 



Gxma{r, 6) e^"^°'^rdrde 
,p, m, a) rdrdd 



For particular characterizations of S'^^{c,p,q,m,a), the resulting patterns of receptors 
can be simulated. These representations can be obtained by means of the computer 
graphics technique of ray tracing [42] using grey tones corresponding to the values of 

Cxmair, 0). 

Results and discussion 

Anderson et al. [7] proposed that LDL receptors in human fibroblasts are inserted at 
random locations and are transported by diffusion toward coated pits, and Goldstein 
et al. [21,25] have studied theoretical aspects of the dynamics of the LDL receptor on 
the cell surface. One of the questions addressed in these models was whether or not 
the random insertion of LDL receptors into the plasma membrane, followed by pure 
diffusion and using the measured diffusion coefficient for LDL receptors, could give a 
rapid enough aggregation of receptors in coated pits to account for the observed rate 
of removal of LDL from the cell surface. Since an experimentally determined lower 
bound for the forward rate constant kdu+ is (2.3 x 10 ± 1.6) x 10"^° crrp'ls and the 
value obtained from their model is 1.9 x 10" crn^/s, the answer to that question is that 
the hypothesis of random insertion and pure diffusion of LDL receptors to coated pits 
is consistent with experimental observations but just barely. Therefore, even though 
the diffusion-limited forward-rate constant k^iu + calculated by these authors is consist- 
ent with available experimental results on the rate of LDL internalization, these do not 
rule out the possibility that the true value of the mean capture time of receptors 
is actually smaller than r^„, the value predicted on the basis of random insertion of 
receptors followed by pure diffusion to coated pits. Indeed, if LDL receptors on human 
fibroblasts are not inserted uniformly into the plasma membrane, but are inserted pref- 
erentially into specialized regions where coated pits form, and if the resulting reduction 
in the time required for LDL receptors to diffuse to coated pits was significant, the 
conclusion that diffusion of LDL receptors to coated pits is the limiting step in the 
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interaction of cell surface LDL receptors with coated pits could be wrong. Also, diffu- 
sion could be aided by a more active transport in order to increase the rate of 
accumulation of receptors in coated pits. In particular, we have studied here the effects 
of a polarized insertion of receptors as well as those related to a transport mechanism 
set by diffusion aided by a retrograde membrane flow. 

Experimentally determined values for coated pit radius a and the outer radius b of 
the reference region D. are respectively a = 0.10 ±0.05 [^m [43] and ^=1.0 ±0.2 nm; 
these were obtained by using equation (5) with the calculated value for the suitable 
density of coated pits p = Q.?>\ ± 0.09/^w^ [16]. In order to typify values of the funda- 
mental ratio A = vl2D (cf Eq. Bl) we will use a low reference value v = Vq for the 
strength v of the retrograde membrane flow. The value Vq = 9.6 x 10" ^ cmlmin repre- 
sents the maximum of a hypothetical local flow that originates when coated pits inva- 
ginate to internalize the trapped ligand-receptor complexes [44]. As a second reference 
value for the strength v of the retrograde flow we may take Vi = 10"* cmlmin. This 
value gives the approximate rate at which objects on the surface of fibroblasts are swept 
backwards [26,14]. Since Vi is approximately equal to IOOvq, this retrograde membrane- 
flow rate would sustain an extremely fast convective transport in comparison with one 
having a strength Vq. In order to characterize a variation range for the diffusion coeffi- 
cient D, we will consider a normal process for the LDL receptor. This is associated with 
a diffusion coefficient value of Dq = (2.7 ± 0.09) x 10"^ cmllmin, which was determined 
by Barak and Webb [45] using fluorescence photo bleaching recovery (FBR) in an 
experimental arrangement involving LDL receptors diffusing on the surface of human 
fibroblast cells. We can also contemplate a fast diffusion process corresponding to 
= 8.1x10"^ cmllmin, this being the average of the values defining the observed 
range of variation for the diffusion coefficient of the dil (3)- LDL receptor complex on 
blebs induced in the cell membrane. And we could also consider an extreme reference 
value of Dext = 6 x 10" cm'Imin, that is D^^t = 7.40Di. 

The phenomenon of capping is characterized by the formation of large patches of pro- 
teins that passively flow away from the regions of membrane exocytosis. Individual pro- 
teins can overcome this convective transport when they have diffusion coefficients (D) of 
about Dgxt [14]. Moreover, a value for D of -De^f/6 is comparable to the largest estimation 
for the diffusion coefficient of the dil (3)- LDL receptor complex on blebs [45]. For the 
parameter values of the LDL system on human fibroblasts, the calculated value for t^„, 
the mean capture time when receptors move by pure diffusion and are inserted uniformly 
over the entire cell surface, is t^,, = 2.92 min (cf Eq. (11)). 

In the present setting, a uniform insertion-rate function can be acquired by, for 
instance taking a = 0, m = bla and SJJ^(w,0) = c, where c is a positive constant and 
S''p{m,0) = 5^^(w,0) = 0. This representation is known here as a c-uniform insertion 
mode (Figure 2 and Eq. A22). A second description of a uniform insertion mode is 
what we call a pq-un\£orm rate function, in which receptors are inserted solely over the 
regions Dp{m,n) and D^(w, tt); this is realized by setting S^{m,n) = S'^{m,n) = c, 
where c is a positive constant (Figure 3 and Eq. A24). Finally, a description that we 
refer to as cpq-un\£or:m insertion is achieved by sorting receptors over the three disjoint 
regions D^im, a), Dpim, a), D.q(m, a) with S'f{m,a) = S''p{m,a) = S''^{m,a) = c, where 
c is a positive constant (Figure 4 and Eq. A26). Consistently, for the case where no 
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convective flow is present {v = 0) and receptors diffuse normally {D = Dq), all these 
representations indistinctively yield txi^^i = i^du- 

For the case in which convective transport is moderate (v = Vq) and receptors diffuse 
normally {D = Dq), that is, the flow strength to diffusion ratio A has a value A = VqI2Dq, 
calculating r^^a as given by equation (B27) for S'^^{c,p,q,m,a) where insertion is uni- 
form over all of D, we obtained txi^^ = 1-003 t^u- Hence, the present model predicts that 
if LDL receptors diffuse normally, a retrograde membrane flow having strength com- 
parable to Vq would have an insignificant effect on t^„. But for that same convection- 
to-diffusion ratio, we found that a locally uniform insertion mode can prompt greater 
receptor trapping rates than those linked to uniform insertion all over D. But now sup- 
pose that receptor insertion takes place over the regions O^-im, a), flp(m, a) and 0.q{m, a), 
at different constant rates S''^{m,a) = c,S''J^{m,a) = p and S''^{m,a) = q; this results 
in what we term here a c/7g'-locally uniform insertion mode (Figure 5 and Eq. A28). 
Then, if for instance we choose m = 1.5, a = n/2, Sdm, a) = 0.10, 5p{m, a) = 0.30 and 
Sq{m, a) = 0.60, the model yields ixma = 0-49rrf„, which amounts to an important reduction 
of Trf„. For a = n we obtain what we label here as a /7^-locally uniform insertion form 
(Figure 6 and Eq. 30) denoted by the symbol S'^^{0,p, q, m, n). Receptors are inserted in a 
radially symmetric manner over each of the regions Oq(m, n) and Qp(m, n), but through 
different constant rates S''p{m, a) = p and S''^{m, a) = q. For m = 2.0, Sp[m, n) = 0.20 and 
8q{m,n) = Q.?>Q, the model yields r^ma = 0.4lTrf„. Now, since for a pq-locdMy uniform 
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insertion form we have dp{m, n) + Sq{m, n) = 1, then letting dp{m, n) approach zero, 
receptor insertion will be gradually accommodated within the region flq{m, a) so that 
eventually, when Sp{m, n) vanishes, we will obtain what we identify as a g'-plaque form 
insertion mode denoted by means of the symbol S''^(0,0, q, m, n) (Figure 7 and 
Eq. A32). This is actually the plaque-form insertion mechanism envisioned by Wofsy 
et al. [40] for modeling preferential insertion as conceived by Robeneck and Hesz 
[39]. For X = VqIDq, and m = 2.0, the mode 5*^^(0, 0, q, m, n) yields Txma = 0.26rrf„, again 
a noticeable reduction in rj„. This mode is radially symmetric and can be already con- 
sidered as a form of polarized insertion. Similarly, if we initially chose a locally uni- 
form mode S'^^(0,p,q,m,n), then letting 8q{m,n) approach zero we will force 
receptors to be mainly sorted over the flpirn, n), region and in due course when Sq{m, n) 
vanishes will produce what we call a jj-peripheral insertion mode symbolized by 
means of S'^^(0,p,0, m, n) (Figure 8 and Eq. A34). This insertion rate function is a 
peripheral form of polarized insertion; receptors are inserted in an annulus contiguous 
to the outer boundary of £1. For A = Vq/Dq, peripheral insertion yields Tx^,^ ^ 1-003 
with the lower bound of 1.003 rj„ attained in the limiting case when m approaches 
one. Hence, if enhancement of LDL receptor trapping rate is required, a peripheral 
insertion mode for m > 1 turns out to be an inefficient mechanism. For instance taking 
A = VqIIDq, and m = 9.7 we obtain Txma = l-15rrf„. 

We can deal with three different non-radially symmetric polarized insertion forms. 
One is obtained if we let S'f{m,a) = 0 and S''q{m,a) = 0, then all recycling receptors 
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Figure 4 The cpq-uniform insertion mode. A device symbolized liere using S (c, c, c, m, a) and is obtained 
by setting (m, a) = Sp{m, a) = S^im, a) = c, witli c a positive constant. 



will be sorted over the region Op{m, a) and we will be dealing with a jj-polarized recep- 
tor insertion device denoted by means m, a) (Figure 9 and Eq. A36). 
Furthermore, for a jj-polarized mode, by reducing the area of the insertion region flp{m, a), 
we could accommodate all recycled receptors in a favorable position relative to 
the cross section of the coated pit perpendicular to flow streamlines. Then, we may 
conjecture that this will allow the effect of convection to be maximal for receptor 
trapping rate enhancement. But we found that this device could at best yield a nar- 
row variation range of r^^a relative to rj„ and not a relevant reduction in this value. 
In fact, for a /^-polarized insertion mode with 0 < a < nl6 and 1.1 < m < bla, we found 
l.Olr^M < Tj^ma ^ l-16r^M. Now, if we assume that 0<a<n, and also that S'f{m, a) = 0 
hold, then receptors will be sorted over the regions flp{m, a) and Oq{m, a), resulting in 
another form of a non-radially symmetric polarized insertion mode. This is regarded here 
as a /7^- polarized receptor insertion form and denoted through S^^{0,p, q, m, a) (Figure 10 
and Eq. A38). In particular, for m = 2.2, a = nl2, Sp{m, a) = 0.20 and Sq{m, nil) = 0.80, we 
calculated Txma = 0.45 rj^. That is, a jcg'-polarized receptor insertion form could lead to a 
significant reduction of r^j,. But if we now set S''f{m,a) = 0 and S''^{m,a) = 0, then all 
recycling receptors will be delivered over the region flq{m, a), and we will have a third 
form of non-radially symmetric polarized insertion mode. This is named ^-polarized 
insertion and is symbolized by means of S''''(0, 0, q, m,a) (Figure 11 and Eq. A40). 
Additionally, this arrangement along with a slow convective transport {v = Vq) and a 
normal diffusion process {D = Dq) could potentially induce a major reduction in t^u- 
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In fact, if we set m = 2.0 and a = n/6 we obtain Txma = 0.26 r^„, which coincides with the 
plaque form insertion mode mentioned above. Hence, for the case A = vq/Dq either 
polarized insertion over the regions fip{m, a) and f2q{m, a) or over the region a), 
both seem to provide an efficient mechanism for the reduction of t^u- 

If we now assume that A = ViI2Dq, that is, a relatively fast convective transport (v = Vi) 
and a normal diffusion process {D = Dq), and if we again choose the ^-polarized inser- 
tion mode S'^^(0, 0, 1, 2, 77/6 ) described above, then we will get r;i,„Q. = 0.83 tj„. This 
result endorses the view that whenever convection becomes a relatively more active 
transport than diffusion, the process will fail to induce an effective randomization in 
the surface distribution of receptors [14]. The relative dominance of convection over 
diffusion implies that those receptors favorably distributed in an area determined by 
the cross section of the coated pit perpendicular to flow streamlines could be trapped, 
while the others would be swept away by the flow. Therefore, whenever convection is 
relatively fast (v = Vi) we may assume that increasing D, that is, taking D > Dq could 
induce relatively greater aggregation rates in coated pits. But since we have -Dq = (2-7 ± 
0.09) X 10"^ cm2lmin, then the largest feasible experimentally determined value for Dq 
is only D = 2.7910"^ cm2/min, and this along with the 5*^^(0, 0, 1, 2, 7t/6 ) mode and a 
fast convective transport will produce T;ima = 0.75 r^a, which only amounts to a moder- 
ate reduction of the value Txma = 0.83 rj„ calculated above. Further, if convection is 
relatively fast (v = Vi) and we still have a ^-polarized insertion mode S'^{0, 0, 1, 2, n/6 ), 
then unless we have D = 1.74Do, the calculated value r^^o: will be significantiy reduced: 
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Figure 6 The pq-locally uniform insertion mode. A receptor insertion meclianism wliicli is denoted by 
means of S'^(0,p, q, m,n). It is radially symmetric over each one of the regions Qp and fiq. This form is 
associated with the case Sf(m,n) = 0, S'^{m,n) = p, and S'^(m,n) = q, with p and q positive constants. 



for that arrangement we obtained T^ma = 0.24 Tdu- Therefore, whenever convection is 
relatively fast {v = vi), even for an efficient insertion mode such as the g'-polarized form 
S''^(0, 0, 1, 2, n/6), larger values than D = I.IADq will be required for a noticeable reduction 
in Trf^. Indeed, for a fast retrograde flow (v = Vi) and a fast diffusion process {D = Di), we 
found that S''^(0, 0, 1, 2, n/6 ) can produce a value of r^ma = 0.01rj„, which is a drastic re- 
duction of Tdu- If we assume, moreover that receptors are returned to the cell surface 
through a /^-polarized insertion mode with 0 <a <n and with 9.5<m< bla, then even if 
convection is extremely fast (v = 300vi), letting D = Dg^i, we get 0.19rrf„ < Txma ^ 0.65t^„. 
Compared with the results obtained above for the same insertion mode, this last inequal- 
ity highlights the importance of diffusion in the reduction of the mean capture time of 
LDL receptors by coated pits. However, comparing with the expected variation range for 
Dq, D values beyond D = 2.7910" * cm2/min seem to be extremely high, and likely not 
achievable in the experimental system analyzed. Thus, before a faster diffusion process 
can be raised as an agent that can induce an effective randomization of the surface distri- 
bution of LDL receptors in the presence of a fast retrograde membrane flow, we assume 
instead that there must be another factor that could be adduced as a feasible device for 
promoting higher aggregation rates of LDL receptors in coated pits. And so, in accordance 
with the results presented here, we can suggest special forms of the insertion rate function 
S^^icip, q, m, a) as mechanisms that potentially could boost the rate of removal of recep- 
tors by coated pits. Thus under the assumption of a fast convective transport, v = v^, and 
a normal diffusion process, if we allow returning receptors to be sorted by means of a 
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g'-polarized insertion form by setting m = 1.3 and a = jil2, we would get ri^„ = 0.23r^„, 
implying a noticeable reduction in rj„. This means that whenever convection is fast and 
receptors diffuse normally, increasing the receptor-trapping rate could depend decidedly 
on polarized receptor insertion at preferential regions near coated pits. 

For the insertion-transport mechanism for LDL receptors under consideration, we 
also address here the theoretical exploration of their consequent display on the cell 
surface. A capping-like surface display should manifest as a graduated concentration of 
receptors in the direction of flow streamlines. Figure 12 displays the surface aggregation 
patterns of receptors that result when a polarized reinsertion mode is considered. The 
combination of polarized insertion, fast convective transport of rate v = Vi and a rela- 
tively slow diffusion process such as one associated to Do can be observed to produce a 
marked gradient in the distribution of unbound receptors (Figure 12a and b). This is 
consistent with the conceptual model for the capping phenomenon presented by 
Bretscher [14]. Further, and in concurrence with that author, the present model 
explains this gradient as being induced by a relative dominance of convection over 
diffusion; that is, an arrangement which precludes an effective randomization of the 
surface-receptor distribution. Thus, in the presence of a retrograde membrane flow 
with typical strength v = Vi, even single LDL particles having a diffusion coefficient 
value D = Do — which is considered normal — would produce a capping-like cell surface 
distribution; thus precluding both a uniform distribution and the display of clusters of 
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Figure 8 The p-peripheral insertion mode. An insertion form denoted here by means of S (0, p, 0, m, n). 

It is a radially symmetric and polarized and linl<ed to the case Sf{m, n) = 0, Sp {m, n) = p, and SJ^(m, 77) = 0, 
witli p a positive constant. 



unbound receptors on the surface. Even in the presence of a fast convective transport, 
a suitably fast diffusion process can reverse the capping effects (Figures 12c and d). 

Our results are mainly consistent with the capping phenomenon, and suggest that 
this is indeed concomitant with the idea of a retrograde membrane flow. Capping is 
described as the rearward transport across the leading lamella of various materials used 
to mark the cell surface [46,14]. Bretscher [14] asserts that during capping, those recep- 
tors that are excluded by coated pits, and therefore do not recycle, may be continuously 
swept by a retrograde membrane flow. This author also states that whether such a non- 
circulating protein does get swept backwards or not depends on how fast it can diffuse 
by Brownian motion and how fast the flow is; and moreover that for individual proteins 
to overcome the sweeping effect of the flow, they must have diffusion coefficients (D) 
of about Dgxt [14]. This assertion was corroborated by our results, since our model pre- 
dicts that for D values of the order of D^^t the capping effect can be reversed. Neverthe- 
less, a value of Dg^t lies far beyond the maximum experimentally feasible value for Dq. 
And even though Bretscher [14] claims that FBR often underestimates the true value of 
the diffusion coefficients of pertinent surface receptors, a value Dg^j for D in the LDL 
experimental system seems unfeasible. Our results suggest that this is where the rela- 
tive importance of non-uniform insertion becomes evident as a factor for boosting LDL 
receptor trapping. Indeed, even when convection is fast and receptors diffuse normally, 
our findings indicate that special forms of polarized receptor insertion can be arranged 
so as to produce a substantial reduction in r^^. This was more obvious in the case of a 
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g'-polarized mode, which could be raised as a feasible preferential-insertion mechanism 
for enhancing the LDL receptor trapping rate. 

The problem of the surface distribution of unbound receptors is an important ques- 
tion pertaining to the LDL experimental system. For example, Wofsy et al. [40] argued 
that in the Robenek and Hesz [39] experiments, the LDL-gold particles were highly 
multivalent and thus may have been bound more efficiently to aggregate than single 
receptors, suggesting that the aggregation of newly inserted LDL receptors in regions 
around coated pits was still unproven. But beyond the caveats of Wofsy et al. [40] the 
results of Gross and Webb [47] were considered to provide a quantitative basis in 
support of a surface distribution of both LDL particles and their receptors in different- 
sized clusters. And again, the existence of these clusters was questioned by Sanan et al. 
[48] who reported experimental results that detected a dispersed or scattered popula- 
tion of LDL receptors, in addition to and clearly distinct from, clusters that formed 
soon after their insertion during recycling. In view of these disparate reports, Robenek 
et al. [49] performed new experiments on fibroblast and hepatocyte plasma membranes 
which confirmed their original claim that the initial insertion and display of LDL recep- 
tors in fibroblasts occurs predominantly as plaques. But experiments using fluores- 
cently labeled LDL and influenza virus particles bound to the surface of human 
fibroblasts imaged with a cooled-scan CCD camera attached to a fluorescence micro- 
scope were found to yield a reasonably accurate measure of the proportion of single 
particles, but large errors were encountered in the proportion of large cluster sizes 
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Figure 10 The pq-polarized insertion mode. This insertion meclianism is denoted liere througli 5 (Q,p, 
q,m,a). It is polarized and non-radially symmetric and linl<ed to tlie conditions Sf{m,a) = 0, Sp{m,a) = p, 
and 5'^(m, a) = q, with p and q as positive constant. Even for )\ = v^/2Da, which produce capping-lil<e 
effects, certain characterization of this mode can give substantial reductions in j^^. 



[50]. The high proportion of single particles found in these data provided evidence 
against the clusters claimed by Robenek et al. [49] and support the random insertion 
model. It is still a matter of disagreement among different research groups whether 
newly inserted LDL receptors are dispersed in the plasma membrane of cultured fibro- 
blasts prior to entrapment in coated pits or whether they remain after reinsertion as 
groups of clusters. Our findings reveal that in the presence of a fast retrograde mem- 
brane flow, the expected surface display of unbound receptors must be consistent with 
the capping phenomenon; thus, the formation of steady-state surface clumps of LDL 
receptors is unachievable in the current setting (Figure 12). 

Traditionally the existence of a retrograde flow has been a source of controversy 
among researchers, some of whom disclaim this paradigm (e.g. [51,52,53]) while 
others convincingly address it (e.g. [14,54,55]). In particular, our results provide a 
theoretical support for the claim of Bretscher [14] which suggests that, if a retrograde 
flow does exist, the diffusion-convection ratio could be a fundamental determinant of 
the surface display of LDL receptors that are not bound to coated pits. And even 
though biochemical or biomechanical factors have been invoked in order to put for- 
ward other paradigms for the capping phenomenon (e.g. [49,56-58]), our findings 
support the views of Bretscher [59] and Ishihara et al. [60] on the adequacy of the 
retrograde-flow model. Hence, the involvement of diffusion, polarized reinsertion and 
a retrograde membrane flow in the formalization of the receptor-mediated endocytic 
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Figure 11 The q-polarized insertion mode. A receptor insertion paradigm denoted liere by means of 
S'®(0, 0, q, m, a] is linked to the case Sf(m, a) = 0, Sp (m, a) = 0, and S!^(m, a) = q, with q a positive constant. 
This mode is non-radially symmetric and polarized. Setting = n , the q-polarized insertion mode gives the 
paradigm of Wofsy et al. [35] for insertion in plaques. Even when A = v-,/2Do holds so that the capping 
phenomenon could create a graduated distribution of unbound receptors in the direction of flow streamlines, 
this mode seems to provide a highly efficient form of receptor insertion by making dramatic reductions on 1^^. 



cycle may provide a trustworthy paradigm for the interpretation of observations. More- 
over, the fact that some characterizations of the retrograde flow-diffusion ratio could in- 
duce slower receptor trapping rates cannot provide a rebuttal for the retrograde-flow 
model. Indeed, Echavarria-Heras et al. [41] suggested that receptors could be projected to 
the plasma membrane in a preferential mode relative to the flow streamlines. We have 
studied here the consequences of such an arrangement and our results show that cer- 
tainly, even in the presence of a fast retrograde flow, we could consider a ^-polarized form 
as a mechanism for the enhancement of LDL receptor trapping rates. 

Conclusions 

Our results seem to indicate that even in the presence of a fast retrograde membrane 
flow, which could influence capping-like effects, diffusing LDL receptors that are sorted 
through polarized insertion in the form conceptualized here, could spur faster aggrega- 
tion rates in coated pits than those associated with uniform insertion over the entire 
cell membrane. Moreover a g'-polarized form could provide a paradigm for a highly 
efficient form of preferential insertion. The effectiveness of this mode is mainly 
explained by its placement relative to the flow streamlines. Furthermore, if the cross 
sections of the linked insertion region flq{m, a) and those associated with a coated pit 
are comparable, then even though a relative dominance of convection over diffusion 
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Figure 12 Surface aggregation patterns of unbound LDL receptors associated with different values 
of the fundamental ratio X and non-radially symmetric-polarized insertion modes, a) the surface 
pattern made by A = v^/IDq and the q-polarized insertion mode S'^(0, 0, q, m, a], with m = 2.3, 6g{m, a) = 1 and 
a = n/4; b) surface pattern associated with the case X = v-,/2Di and the p-polarized insertion mode S'®(0, p, 0, m, a) 
with a = 77/5, 771 = 9.5, and 6p(m,a) = 1; c) the surface aggregation pattern formed by A = \/,/2D„f 
and q-polarized insertion S'^{0,p, q, m, a) with m= ].2, a = 77/5 and 5^(771, a) = 1; d) the surface receptor 
pattern for X = v-,/2D^),t and p-polarized insertion S'^(0, p, q, m, a) with m = 9.5, a = 77/4 and 6p{m, a) = 1. The 
patterns shown are consistent with the capping phenomena; that is, when convection is fast iv = v-,) and 
diffusion is normal (D = Do) a graduated distribution in the direction of flow streamlines is observed, a) 
and b). But in the presence of a fast convective transport, suitable values of the diffusion coefficient 
(e. g.D = Dext) can reverse the capping effect and induce an effective randomization of the distribution of 
unbound LDL receptors, c) and d). In any event, in the presence of a retrograde membrane flow, no 
receptor surface clusters are formed. 



might occur, the privileged position of the insertion region Oq{m, a) would induce a 
fast receptor-trapping rate. But our results also hint that in spite of the efficiency of this 
device, the influence of a retrograde flow would cancel out the display of the corre- 
sponding surface clusters. Moreover, if new experiments do reveal that the surface 
clusters exist, then it could be necessary to envision a mechanism that explains their 
formation in the presence of a retrograde membrane flow. Meanwhile, for the case 
V = Vi and D = Dq, a ^-polarized insertion mode with m = 1.20 and a = it/ 16 provides a 
convenient paradigm for explaining the coexistence of a retrograde membrane flow and 
receptor aggregations that might be experimentally misinterpreted as surface clusters. 
Indeed, considering the experimental error associated with the measurement of the 
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coated pit radius, a region of insertion of the form with m = 1.20 could actually be 
a region placed over the boundary of a coated pit. And the surface clusters observed 
might be formed when recycling receptors are returned near the periphery of newly- 
formed coated pits so as to enhance the removal of the LDL ligand, but the existence 
of that insertion paradigm would require experimental substantiation. 

Appendix A 

In accordance with Eq. (17), the function S'^^{c,p,q,m,a) inserts S'f{m,a) particles per 
unit area per unit time at (r, 6) in D.. Suppose that the number of particles that S''^{m, a) 
inserts in Q' is L{m, a). Then, we have 

fb /■TT-a 

Sf (m, a) rdrdd = /, {m, a), (Al) 



a J -n+a 
b pjT+a 



/ / Sf{m,a)rdrde = Ip{m,a), (A2) 

./ ma J TT-a 
ma p7T+a 

S'^im, a) rdrdd = Iq{m, a) . (A3) 



a J n-a 



Let's denote by means of /q the total number of particles inserted by S^^{c,p, q, m, a) 
in D.. Assume also, that a proportion Sp(m, a) of the total number of particles inserted 
in 13 is sorted by S'^^{c,p, q, m, a) in D,p{m, a), that a proportion Sq{m, a) of In is placed 
into the region 0.q{m, a), and that a proportion Sdm, a) of In is delivered in D,q{m, a). 
Then, we must have 

Ic {m, a) = Sc{m, a)In, (A4) 
Ip{m, a) = Sp{m, a)In, (A5) 
Iq{m, a) = dq{m, a)In, (A6) 

and the balance equation 

8p{m,a) -\- 8q{m,a) -\- 8c{m,a) = \. (A7) 

From equations (Al) through (A3) and the mean value theorem there are numbers 

(r„0,) e a(w,a), (A8) 

{rp.dp) e Dp(w,a), (A9) 

[rq.dq) &nq{m,a), (AlO) 
such that 

aS';''''{m,a,){b^-a^m^) 

6p{m, a) = , (All) 

5q(m,a)= , (A12) 



8c{m,a) 



aSf':{m,a,){n-a){b^-a^) 



(A13) 
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Therefore, we have to introduce the consistency conditions 



8p{b/a, a) = 0, 


(A14) 


SJm,0) = 0, 


(A15) 


Op[l,n) = 1, 


(A16) 


dq{l,a) = 0, 


(A17) 


Sq{b/a,n) = 1, 


(A18) 


dq{m,()) = 0, 


(A19) 


Sc{m,n) = 0, 


(A20) 


Sc{m,0) = 1. 


(A21) 



A uniform insertion mode is obtained by setting S''"{b/a,0) = c, where c is a positive 



constant. This is symboUcally represented by means of S'^^{c, 0, 0, b/a, 0) and given by 

(c {r,d)en,{b/a,Q) 
S'^c, 0, 0, b/a, 0) = <^ 0 (r, B) e Dpib/a, 0) , (A22) 

{r,e)£nq{b/a,0) 



where 

In 



n{b- 



(A23) 



Also, uniform insertion is associated in the case S''f{m, n) = 0, S'^{m, n) = S'^{in, n) = c. 
This mode is known as /(^-uniform insertion mode and denoted by means of 
S''^(0,p, q, m, n), 

(0 {r,d)sn,{m,n) 
S'^iO, c, c,m,n)= I c (r, 6) e Dp(w, n) . (A24) 

Again, from equations (All) through (A13), 

In 



a{b^-a^) ■ 
Moreover, the insertion rate function 



(A25) 



{c {r,d) e nc{m,a) 
c {r,8) snp{m,a) , (A26) 
c {r,d) snq{m,a) 



which is obtained for the case Sfim^a) = S''^{m,a) = S''^{m,a) = c, also inserts parti- 
cles uniformly all over O. This is known as a c/jg'-uniform insertion mode. From equa- 
tions (All) through (A13) we have 



a{b^-a?-rrfi^ -\- {n-a){b^-a?-^ -\- a{a?-rrp--a^ 



(A27) 
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The form of S''^{c,p,q,m,a) obtained by setting S''f{m,a) = c, S^p{m,a) = p and 
S'^{m,a) = q, where c, p and q are positive constants will be known as a cpq-locally 
uniform insertion mode. That is, from equation (17) we will have 

{c {r,d) e nc{m,a) 
p (r,6l) eDp(w,a) (A28) 
q {r,e) e D^(w,a) 

with 

apip^-a^m^) ^aq{a^np--a^) ^ac{n-a)(b^-a^) ^ (A29) 
In In In 

A second form of locally uniform insertion denoted by means of 5''''(0 , p, q, m, n) and 
known as a /7^-locally uniform insertion mode is obtained for the case S''f{m,n) = 0, 

are 
1 



S'^{m, jt) = p and S''^{m, n) = q, that is, 



(0 {r,d)sn,{m,n) 
S'\0,p, q,m,n) = Ip (r, d) e X2p(w, n) (A30) 
\q (r,0) eX2^(w,7z-) 

where, from equations (All) through (A13) p and q satisfy 

^ + ^^^— ^=1. (A31) 

la In 

The case in which receptors are inserted into two disjoint regions of abruptly con- 
trasting rates is defined as a polarized insertion mode. It will be typified here as a mode 
that sorts the total number of recycled receptors over either the region D.p{m, a), or 
Dq{m, a), or flp(m, a) u Qq{m, a). Moreover, a polarized insertion mode can be obtained 
in five different forms, two of these forms being radially symmetric. 

The insertion rate function known as a ^-plaque form is denoted by means of 
S''^{Q,0,q,m,n) and obtained for S[%m,n) = 0, S'^{m,n) = 0, and S'^{m,n) 
being q a positive constant, that is. 



^^(0, 0, q,m,n) = I 0 (r, 6) e ^2^(w, a) (A32) 
yq {r,d) s nq{m,a) 

is radially symmetric and polarized, and from equations (All) through (A13) we have 
In 



ji{a?-rrp--a?- 



(A33) 



The function S'^O.O, q, m, n) moreover, coincides with the plaque form insertion 
mode proposed by Wofsy et al. [40]. And, by setting 5''^(0, 0, q, m, n) = Sr'^ one obtains 
from equation (17) the decreasing and radially symmetric insertion-rate function 
considered by Echavarria-Heras et al. [41]. 

Further, what we call here a /^-peripheral insertion mode was proposed by Bretscher 
[38]. It can be also considered as a radially symmetric polarized insertion mode. It is 
denoted through S'^^{0,p, 0, m, n) and is obtained by setting S'f{m,n) = 0, S''p{m,n) = p 
and S'^{m, n) = 0, with p a positive constant, that is 
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(0 ir,d)en,{m,n) 
S''\0,p, 0,m,n) = Ip (r, 6) e Dpim, n) , (A34) 



and from equations (All) through (A13) we have 
In 



(A35) 



A first form of a non-radially symmetric polarized insertion mode is known here as 
jtJ-polarized. It is denoted by means m, a) and arises whenever S'f{m,a) = 0, 

S^{m, a) = p and S'^{m, a) = 0, where is a positive constant, that is, 

rO {r,e)&a,{m,a) 
S'%0,p, 0,m,a) = Ip (r, 6) <= Dpim, a) (A36) 



where from equations (All) through (A13) we have 
In 



(A37) 



A second form of a non-radially symmetric polarized receptor insertion mode is 
known here as /?^-polarized. It is denoted by means of S'^^{0,p, q, m, a), and is linked to 
the case S'f{m,a) = 0, S''p{m,a) = p, and S''^{m,a) = q, with p and q positive con- 
stants, that is, 

(0 {r,B) sn,{m,a) 
S''%0,p, q,m,a) = Ip (r, 6) e Dpim, a) (A38) 
[q {r,d) G nq{m,a) 

and from equations (All) through (A13) we have 
Inp{m, a) Inq{m, a) 



a{b^-a?-nf-) a{a?-rrp--a'^ 



1. (A39) 



A third non-radially symmetric and polarized form called ^-polarized is symbolized 
by means of S'^{0,Q,q,m,a). It is linked to the case S'f{m,a)=Q, S'^{m,a) = 0 and 
S'^{m,a) = q, where g' is a positive constant, that is, 

r 0 {r. d) snc{m,a) 
S''%0, 0, q,m,a) = <0 (V, 6) e Dp(w, a) (A40) 
[q {r,d) s nq{m,a) 

where from equations (All) through (A13) we have 

q = , ^\ 2s ■ (A41) 



Appendix B 

In this section we outline the adaptation of the model for receptor mediated endo- 
cytosis presented in Echavarria-Heras et al. [41] for the case of a partitioned recep- 
tor insertion mode. For detailed mathematical proofs the reader is referred to 
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Echavarria-Heras and Leal-Ramirez [61]. In the present settings the reference 

region D, = l^{x,y)\a <\J + >'^< ^| coincides with the Berg-Purcell [23] geometry. 

Equation (5) scales the outer radius h. A retrograde membrane flow will be represented by 

means of the vector field V (x, y. t) = v i with i being the unitary vector pointing in the 
direction of the positive x-axis and v, the strength of the flow in units of cm/s. Since we 
have assumed here a convective diffusion transport and also a partitioned receptor- 
insertion mode (cf Eq. (17)) we must extend the model of equations (20) through (22) in 
order to include these mechanisms. The symbol D will stand for the 2-dimensional diffu- 
sion coefficient, which has units of cnP'ls and will be assumed to be constant. The ratio 

A = — , (Bl) 
2D' ^ ' 

having units cnf ^ is known as the fundamental ratio of the transport process because 
it combines the parameters v and D, which represent the factors that control receptor 
movement on the cell surface. Also, S{m,a,x,y) will denote the rate of particle inser- 
tion, with units: particles/ crn^ - s. At steady state, the corresponding concentration of 
particles C{X,x,y) satisfies the partial differential equation in the annulus D [62], 

-div{DWC) + divide = S{m, a,x,y) (B2) 
Equivalently, 

DV^C-v— + S{m,a,x,y) = 0. (B3) 
ax 

Due to the geometry determined by D, it is convenient to switch to polar coordi- 
nates. Then, if we denote by means of n the normal unit vector pointing out radially 
towards (r, 6) for 0 < 0 < 277, and if / stands for the flux vector, the projection of the flux 
vector over n is given by 

— * riC 

J ■n=D— + vCcosd. (B4) 
ax 

We will consider the partitioned insertion mode of equation (17), that is 

[Sf(m,«) (r,0)e^2,(m,a) 
Sim, a, r,d)= I Sf{m, a) (r, 6) e Dp(w, a) . (B5) 

Because of the nature of the convective vector field Cxmair, 0) is not expected to be radi- 
ally symmetric if the strength of the membrane flow is different from zero. However, 
CxmJj, 6) wUl be symmetric about the x-axis. This property is formally expressed by 

Cxma{r, 0) = Cxma{r, -0). (B6) 

The boundary at r = a remains as absorbing, i.e. 

CAUr,e)\r=a = 0. (B7) 

On the other hand, the reflecting boundary condition given by equation (8) must be 
modified in order to account for the influence of the retrograde flow. We will consider 
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that this boundary condition can be replaced by a flux-vanishing boundary condition at 
r = b. Hence the boundary condition at r = b will be formally represented by means of 
the equation 

J -n d=0. (B8) 

Where dD.b denotes the outer boundary of the annulus D,, since we have assumed a 
steady-state for the receptors. This requirement for a vanishing flux of particles across 
the boundary at r =b can be considered a model for a dynamics in which those 
receptors crossing the boundary at r= b are transported into a contiguous influence 
area determined by the direction of the flow streamlines. 

The further formal treatment can be greatly simplified by introducing a function 
Gimair, 0) defined through, 

Gxma{r,e) = C,ma{r.e)e-^'-'°''. (B9) 

The substitution of Eq. (B9) into Eq. (B3) and the fulfilment of the symmetry and 
boundary conditions for CxmAri 6) show that, Gxinai^i 0) rnust satisfy 

^ Gi,ma(r,o)-A Q7i,ma(r,0) = — 6 , (BlOj 

Gx^aia,d) = 0, (Bll) 

Gxma{r, 6) = Gxmair, -6), (B12) 

Gxmair, 6) periodic in d, (B13) 
U^Gxmair. 6) 



-A COS e Gxynaib: e))e-^'"°"'dd = 0 



(B14) 



Finally, for consistency a limiting radial symmetry condition must be added when the 
Strength of the flow approaches zero. That is, we will require that for fixed D>0, 
Gxmair, 0) satisfies 

Vimx-^oGxma{r, 6) = Qs(r) (B15) 

Where Os(r) is the solution of equation (6) subject to boundary conditions (7) and 
(8). 

Using the boundary conditions (Bll) trough (B14) it is possible to obtain, using 
separation of variables, a solution GxmJj'i S) for equation (BIO). Furthermore GxmJft 0) 
can be written in the form, 

GxUr. B) = YAmair) coskO, (B16) 

/f=0 

where for /c = 0, 1 a<m< bla, 0<a<n and AeR, the function gx^ai'^) satisfies the 
non-homogeneous problem 

L\{u)+X''u = ^ I S''^{c,p,q,m,a)e-^''°'^ cos kBdd, (B17) 
with 
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1 

77 



/c= 1,2,..., 



and 



Z,^(m) = -( -(«') +— m) (B19) 



r \ ( ) ~'~ r 

defined for m e C^([a, ^]) and satisfying boundary conditions 

u{a) = 0, (B20) 

u{b) + Yi,{X)u{b) = 0. (B21) 

Where is the function given by 

C -Xh{Xb)lh{Xb) X*Q 
YkW = { : (B22) 

[ -k/b X = 0 

with being the modified Bessel function of the first class and of order k for k = 
0, 1 and X sR. Notice also that the operator L'^{u) given by equation (B19), as well 
as, the boundary conditions (B20) and (B21) originate when eigenfunctions (p^ir, 9) 
for the operator = - are supposed to have the form (f)x[r, 0) = Rx{r)il/{d). The 
function YkW given by equation (B22) is motivated by the requirement of equation 
(B15) of radial symmetry in the limiting case A — > 0 for D > Q fixed. 
For the operator L\{u) a Green's function H\{r. t) exists. It is given by 

[h{Xa)Kk{Xr)-Kk{Xa)h{Xr)]h{Xt) 

r<t 



"'^''^^ ^ [h {Xa)Kk jXt) - A {la)h (A^)]4 jXr) ^^^^^ 

h{Xa) 

where /<^(-) stands for the modified Bessel function of the second kind and of order r. 
Using formulas (9.6.10) and (9.6.53) in Abranowitz and Stegun [63] we can obtain an 
expression for H\{r, t) depending only on 4(z). This permits the proof of the existence 
of the solution Gxmai^t 9) in form given by equation (B16) and which depends continu- 
ously on the boundary data, namely 

C,ma{r. 6) = Y^Lair)<^''°'' COS kd (B24) 

with ^x,y^^{r) expressed in terms of the Green's function H\{r, t) through the expression 

gLair) = ^j'(^-tHl{r,t) jy'{c,p,q,m,a)e-'"^''^' cos kOdd ^t, (B25) 

and ai, given by equation (B18) [61]. 

The symbol kxma+ denotes the pertinent forward rate constant, which can be ob- 
tained by calculating the number of receptors hitting the trap in unit time and dividing 
it by the average concentration (cf. Eq. (1)). We have nevertheless assumed the 
existence of a steady-state concentration Cxmai^i Q) in the reference annulus D, for all \. 
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So the number of receptors hitting the trap in unit time coincides with the number of 
particles that are inserted into the reference annulus in unit time. Hence, we have 
equivalently 

I I S''^{c,p,q,m,a) rdrdd 

Ja J-n (^^2e) 



<Xma+ 



^ / Cxma{r,d)rdrdd 

J a J-n 



and by virtue of equations (1) and (3) we must have for Txma the associated mean 
capture time, 



^Amct — „'l (B27) 







[ gLairy^"' " COS kdrdrdd 

-n 


ima — 

i 




S^^{c^p, q, w, a)rdrdO 

-IT 



For the case of a uniform receptor insertion (i.e. the case m = bla and a = 0 in equation 
(B5)), we will have S'^^{c,p, q, m, a) = S, a constant ior a<r<b. Denoting by ^^^(r) the cor- 
responding form oig\^^{r) and then using the result, 

j e^'™' " cos kOdd = 2TT-{-\fh{Xt) (B28) 

(cf Abramowitz and Stegun [63]) equation (B25) yields 

g\s = '"^^1"'^'"' ftH\[r, t)h{Xt)dt. (B29) 



D 

Meanwhile, the corresponding mean capture time r^^ becomes 

Er=o^'^(-l)'^' f tH\{r,t)h{Xt)dt\h{Xr)dr 
Txs = — ^ , (B30) 

where 

C 2 if k=0 

9k ={ 

[ 4 if k>l 

It can be also shown [61] that for D >0 and fixed, we have 

m, Txs = ; B31 

that is, if receptor insertion is uniform all over D, then whenever A approaches zero, 
Txs approaches the value t^u of equation (12), calculated by Berg and Purcell [23]. 
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